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SUMMARY 


This report documents a variety of related scramjet engine topics in 
which primary emphasis is on simplicity and conceptual clarity. Thus, the 
flow is assumed to be one dimensional, the gas is thermally and calorically 
perfect, and the study focuses on low hypersonic Mach numbers. The first 
technical section evaluates the thrust and lift of an exposed half nozzle, 
which is used on the aero-space plane, as well as one that is fully confined. 
A rough estimate is provided of the drag of an aero-space plane. Background 
material dealing with thermal effects and shock waves is discussed in the next 
section. The following section then presents a parametric scramjet model, 
based on the influence coefficient method, that evaluates the dominant 
scramjet processes. The independent parameters are the ratio of specific 
heats, a nondimens ional heat addition parameter, and four Mach numbers. The 
total thrust generated by the combustor and nozzle is shown to be independent 
of the heat release distribution and the combustor exit Mach number, providing 
thermal choking is avoided. An operating condition for the combustor is found 
that maximizes the thrust. An alternative condition is explored when this 
optimum is no longer realistic. This condition provides a favorable pressure 
gradient and a reasonable area ratio for the combustor. The next section 
provides parametric results based on the model. One significant finding is 
the sensitivity of the thrust to the value of the ratio of specific heats for 
the air upstream of the combustor. The final section summarizes and discusses 
the analysis. 


INTRODUCTION 

In recent years, there has been a renewed interest in scramjet engines. 
These engines operate at hypersonic vehicle speeds, and the flow internal to 
the engine is supersonic. Above a free-stream Mach number of about 5, a 
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ramjet engine is no longer viable. Only a scramjet or rocket engine can 
provide a significant thrust in this flight regime. For powered flight in the 
atmosphere over a long distance, only a scramjet engine is suitable, since the 
ambient air replaces the rocket* s stored oxidizer. 

Conceptually, a scramjet engine is simple, although the problems 
associated with it are severe and some of them are still unresolved. For 
instance, one of these problems is the difficulty of performing appropriate 
tests at hypersonic flight Mach numbers. Because of this difficulty, there 
has been an effort to computationally model the engine and its inlet (see 
references 1 and 2). Over the past 30 years, there have also been a number of 
system-oriented models (see references 3-7). However, it is often difficult 
to assess the validity of these models. They frequently utilize correlation 
formulas and numerous efficiency or fudge factors. Despite these factors, it 
is not evident that the analysis adheres to the governing conservation laws of 
fluid dynamics. 

In this report, a variety of issues are discussed. In each instance, 
emphasis is on simplicity and conceptual clarity. A primary objective of this 
effort is tutorial where we place the scramjet engine within a simple 
theoretical framework. Nevertheless, we are not aware of a similar document 
in the journal literature or among those reports in our possession. 
Consequently, certain aspects of the analysis may be new. 

To lay the groundwork for the subsequent discussion, we evaluate the 
thrust and lift of an exposed half nozzle and one that is fully confined. The 
next section then concludes with a rough estimate of the drag of an aero-space 
plane. This estimate is needed to calibrate subsequent thrust estimates. The 
following section presents background material, which is useful in the 
subsequent sections. Thermal considerations and shock waves are the principal 
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topics discussed. A parametric scramjet model, based on the influence 
coefficient method, is then provided. Although elementary, the model 
evaluates some of the dominant scramjet processes to first order. Results 
based on the model are presented in an extensive series of tables, while the 
final section summarizes our findings. 


SYMBOLS 

a speed of sound 

A cross-sectional area 



I 

sp 

L 


fixed parameters in equation (47) 
drag coefficient 
skin-friction coefficient 
drag 

hydraulic diameter 

unit vectors along the x and y axis 
impulse function 
specific impulse 
lift 


m mass flow rate per unit depth 

m^ fuel mass flow rate 

M Mach number 


n unit normal vector 


p pressure 

q heat addition per unit mass 

Q nondimens ional heat addition 

Q r maximum heat addition in Rayleigh flow 

R gas constant 

S surface area 
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pi an form area 
temperature 
flow speed 
velocity 

Cartesian coordinates 
thrust 

defined by equation (53) 
ratio of specific heats 
Mach angle 
density 

equivalence ratio 
subscripts and superscripts 

B end of kernel region on lower wall (see figure 3) 

c confined nozzle 

d downstream of a normal shock wave 

e exposed half nozzle 

f nozzle exit 

i inlet 

is isentropic 

j lower edge of jet 

mx maximum 

u upstream of a normal shock wave 

w wall 

wJ? lower wall 

wu upper wall 

0 stagnation 

1 combustor inlet 


S 

P 

T 

V 

-► 

V 

x,y 

r 

« 

V 
P 
P 
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1-2 combustor 

1-3 combustor plus nozzle 

2 combustor exit or nozzle inlet 

3 nozzle exit 

oo free stream 

( )* M = 1 condition 

( ) dimensional value 


FORCE DISCUSSION 

Figure 1 shows a conventional sketch of an aero-space plane in hypersonic 
flight. The afterbody consists of what we call an exposed half nozzle. 
Internal to the scramjet engine, a high pressure and temperature gaseous state 
is achieved. This state is a result of some compression by the forebody 
itself and its bow shock wave, by the scramjet’s inlet, and by air/fuel 
combustion downstream of the inlet. Generally, some expansion occurs inside 
the engine before the inlet to the exposed half nozzle is reached, where a 
further expansion occurs. Both expansions, internal and external, generate 
thrust. Later, we evaluate the total thrust. These later sections show that 
internally generated thrust is quite important. 

Our discussion in this section is confined to the exposed half nozzle, 
which produces both lift and thrust. Lift production, on the other hand, does 
not occur in the internal nozzle, where the flow field is approximately 
symmetric. If the vehicle were a rocket, the thrust nozzle would be fully 
enclosed. An exposed half nozzle is used to reduce the vehicle* s weight, 
size, and heat transfer load. We anticipate that the lift may be of secondary 
importance relative to the lift generated by the forebody. 

Our objective in this section is to derive simple formulas for the lift 
and thrust developed by the exposed nozzle. For the subsequent analysis, we 
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need similar relations for an enclosed bell— shaped thrust nozzle. Our 
starting point utilizes the integral relations developed in references 8 and 9 
for waverider vehicles. This approach is most conveniently found in section 
20.3 of reference 10. (In the subsequent discussion, we shall often have 
recourse to this reference.) 

The equations for mass and momentum of a steady, inviscid flow can be 
written as 



-f ~ 

V*nds = 0 


a) 


f [pV(V-n) + pn]ds = 0 
J c 


( 2 ) 


where p, V, and p are the density, velocity, and pressure, respectively. The 
integrals are over a simple closed surface S, which has an outward unit normal 

/V 

vector n. Aside from the previously mentioned assumptions, we also assume a 
two-dimensional flow of unit depth. 

Confined Nozzle 

Figure 2 shows a typical divergent nozzle with uniform inlet and outlet 
flows. The inlet and exit Mach numbers, and M^, are arbitrary, although 

for the subsequent comparison with an exposed nozzle we assume 

1 $ M i < M f (3) 

The inlet and exit areas are y^ and y^, as shown in the figure. The enclosing 
surface consists of three sections, i.e., S^, S^., and S^, where the upper and 
lower wall surfaces are symmetrical. On these surfaces, we have: 
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A A 


V. = V.e , n = 
1 ix 


on S . 


(4a) 


A A A 


V „ - V«e , n = e 
f f x* x 


on S, 


(4b) 


Vn = 0 


on S 


w 


(4c) 


With these relations, equations (1) and (2) become 


m = (pVy) i = (pVy) f 


(5) 


(pv 2y ) , e x - (py)i e x + J s P nds + (p y 2 y) f e x + (py)^ = o 


( 6 ) 


where m is the mass flow rate per unit depth, and e and e are unit Cartesian 

x y 

vectors along their respective coordinates. 

4 

As indicated in figure 2, we assume the upstream velocity and nozzle 


exit velocity are parallel. The lift L and thrust relative to V^, 

then defined by 


are 


L = e 


f pnds 

J o 


Sw 


(7) 


c/= - e • I pnds 
J Sw 


( 8 ) 


respectively. With the use of equations (5) and (6), we readily obtain 


L = 0 

c 


( 9 ) 


^ = (py + mV) f - (py + mV^ = F f - 


( 10 ) 
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where a c subscript denotes a confined nozzle, and the impulse function F is 
defined as 

F = pA + pV 2 A (11) 

As expected, the thrust ¥ is simply the difference between the inlet and 
exit values of the impulse function. 

It is worth noting that equation (10) is not directly applicable to a 
conventional thrust rocket (see reference 11). Our analysis, therefore, is 
not appropriate for this propulsion device. 

Exposed Half Nozzle 

Figure 3 shows a schematic of an exposed nozzle. The inlet and exit 
flows are again uniform and parallel with the same inlet and exit areas as for 
the confined nozzle. The flow is confined by an upper wall and a lower 

wall that extends from the origin to point B on the x-axis. The shape of 

the upper wall can be based on the theory for a two-dimensional minimum length 
nozzle with a planar inlet surface (see reference 10), although this specific 
configuration is not necessary for our purposes. In any case, the flow field 
should be viewed as the upper half of that for a confined symmetric nozzle. 
The flow is uniform downstream of the BC Mach line, which is straight. It is 
simpler, however, to use for the enclosing surface the sum of five surfaces 


s.+s +S^ + S. + S. 
i wu f j wJ? 


where S . is along the exposed lower edge of the jet and 
0 


V^. For this analysis, inequalities (3) are required. 

On S., S„, and S , equations (4) still hold, 

l f wu 

surfaces, we have 


^ is perpendicular to 


On the remaining two 
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V- n = 0, n = 


e on S A and S . 

y w? j 


( 12 ) 


Equations (1) and (2) now yield equations (5) and 


A J» A A r r X B 

- F i e x + J_ pnds + ¥x - Un PdX + P f (X ^ VJ 

‘Vu u 


e =0 

y 


(13) 


The exposed nozzle, denoted by an e subscript, provides a lift given by 


L = e 
e y 


IT pnds] - f 8 pdx 

U Swu J J 0 


where the rightmost integral is along S^. By multiplying equation (13) with 

A 

*e , we obtain 

tr * 


L e = P f (x f ' X B } 


(14) 


Hence, the nozzle has a positive lift. Its magnitude is quite small, compared 
to the forebody’s lift, because p^ is generally small. Additionally, the 
nozzle contributes a pitching moment about the nose of the vehicle, which 
should not be neglected. 

In view of the parallel flow in both the inlet and exit planes, it may 
seem strange that the nozzle has a non-zero lift. However, the same result is 


obtained by summing the e^ components of the forces on the fluid inside the 
closed surface S. By Newton’s second law, these sum to zero. Basically, our 
result stems from the asymmetry of an exposed half nozzle. 

Equation (14) can be put into a more convenient form by observing that 
the B C Mach line is at a Mach angle where 


yf 

x f - x B 


tan |i f = (mJ - l) 


- 1/2 


(15) 
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The area ratio, y^/y*, where by definition M* = 1, is provided by the 
isentropic relation 


h 

y* 


t 2 v(Y+D/[2(T-D] 

lT + li 




(Y+1)/[2(Y-1)] 


(16) 


We have now assumed a perfect gas and a constant value for the ratio of 
specific heats, y. The pressure p^ is similarly given by 


P* = P r 


(i ♦ ^ «?) 


-T(Y-l) 


(17) 


where a zero subscript denotes a stagnation quantity. With the aid of 
equations (15)-(17), our final expression for the lift is 


L 

£ _ j 

( 2 ) 

(Y+1)/[2(y~1)] | 

("1 - !) 1/2 

* 1 
P 0 y 

It + iJ 

M f | 

(l ♦ * «?) 

I 1 " 


(18) 


The thrust of the exposed nozzle is given by 



A 

This time, we multiply equation (13) by -e^, to obtain 



(19) 


Hence, we have the result that the confined and exposed nozzles have the same 
thrust. This conclusion is useful in the subsequent modeling where, for 
simplicity, a confined nozzle is used instead of the actual exposed nozzle 
shown in figure 1. 
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Drag Estimate 


We need a rough estimate of the drag of an aero-space plane in order to 
know if our subsequent thrust estimates are sufficient to overcome the drag 
and provide acceleration. As usual, the drag coefficient is defined as 

C = 2P = 2D 

° (pv 2 )^ S p YPo^Sp 


/v 

where D is the dimensional drag, p and M are the free stream pressure and 

03 00 

Mach number, and is the projected planform area. 

For the later scramjet analysis, it is more convenient to introduce a 
different nondimens ional drag given by 


D = 


A 



PeoAl 


( 20 ) 


where A. is the projected combustor inlet area on a plane perpendicular to V . 

X oo 

The actual combustor inlet area would be larger if it is canted relative to 
V . We now have 

ao 


D = 9 n 2 I 2 c n 

Z ao Aj D 


For our estimate, we utilize 

Y = 1-4 , M = 7 , = 10~ 2 , C = 6 x 10 

00 Op II 


to obtain 


( 21 ) 


D = 206 (22) 

The Cp estimate in equations (21) is significantly higher than would be 
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estimated for the forebody drag, even if a turbulent boundary layer is 
assumed. However, our estimate also includes internal engine viscous drag 
and, more importantly, the drag associated with compressing the engine air 
flow from its upstream condition to that at the inlet of the combustor. While 
the above D value is a crude estimate, it nevertheless is in the range of drag 
values that can be extracted from reference 7. Due to the p A, 

a a ± 

normalization, D should not rapidly vary with flight altitude. It will vary 
with M , since the drag due to engine air compression depends on M . This 

CD CO 

compression, however, also strongly depends on the Mach number at the inlet of 
the combustor. As a consequence, D is not proportional to 

CD. 

BACKGROUND DISCUSSION 
Overview 

We first outline the major assumptions and provide a brief overview of 
the flow field. The flow is assumed to pass through one, or more, weak 
oblique shock waves before it reaches the inlet of the combustor. Some of 
this compression may be isentropic; the amount depends on the vehicle’s 
forebody and engine inlet design. For purposes of simplicity, we assume M 

GO 

and p isentropically change to M.. and p, at the inlet of the combustor. Of 
°° II 

course, this assumption is optimistic; it should be altered when there is a 
known configuration upstream of the combustor. Although optimistic, the 
assumption is not grossly in error as long as the shock waves are weak. 

A schematic of the engine is shown in figure 4. It consists of just two 
units, a combustor and a confined, isentropic nozzle. The flow at each of the 
four stations is assumed to be uniform as required by our quasi— one- 
dimensional approach. As the figure indicates, the cross-sectional area 
derivative, dA/dx, is typically discontinuous at the inlet and exit of the 
combustor. In addition, the flow at station 1 is considered to be sonic or 
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supersonic. The gas is air, taken as thermally and calorically perfect, with 
a constant value for the ratio of specific heats y. 

The combustor, located between stations 1 and 2, has a variable cross- 
sectional area with inviscid flow and with heat addition, which stems from an 
air- fuel combustion process. We ignore the small additional fuel mass flow 
rate, the fuel-air mixing process, and any changes in average specific heats 
or molecular weight. Due to the heat addition, we anticipate that the cross- 
sectional area ratio Ag/A^ will exceed unity. Thermal choking is avoided by 
requiring that Mg 2 1, where 

Between stations 2 and 3, the flow isentropically expands in a nozzle 
from Mg to 


Mg. Later, we use the pressure ratio 


Mg is to be prescribed. 


£3 _ £§ p 2 P 1 _ 


Y/(Y-0 


,Y/(Y-1) 

p 2 

p « " P 2 P 1 P co " 

1 + V "1 



P 1 


(23) 


where Pg/p^ is evaluated in the analysis section and 


£3 

p 2 






r/Or-i) 


(24) 


We also need the nozzle area ratio 


*2 " M 3 


f 1 + ¥ "f 




(Y + 1)/[2(Y~1) 1 


(25) 


A Brayton cycle, which models a jet engine, consists of an isentropic 
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compression, followed by constant pressure heat addition, which is then 
followed by an isentropic expansion. Our engine differs from this cycle in 
that the heat addition need not occur at a constant pressure. In a 
conventional Brayton cycle, where the heat addition is subsonic, the constant 
pressure assumption is realistic. However, for heat addition in a supersonic 
flow, this assumption is no longer warranted. 

Our objective is to determine the thrust ^^-2 P rov ided by the combustor 
and the thrust of the overall system. A second objective is to evaluate 

various pressure, temperature, and area ratios. In view of the parametric 
nature of the study, care must be taken to avoid unrealistic parameter values. 
These ratios assist in deciding this issue. 

Combustion 

We presume a steady combustion process can occur at supersonic speeds. 
Early scramjet studies were unsure on this point. The issue, however, is now 
mute in view of the success of the cw supersonic chemical laser. In both 
devices, fuel and oxidizer streams mix at supersonic speeds. After mixing, a 
steady combustion process occurs providing the local static temperature is 
sufficient to overcome the activation energy of the exothermic reactions. In 
line with our one-dimensional assumption, we henceforth ignore the effects of 
mixing. While this process is important, it is outside the scope of this 
study. 

Because the Hg fuel is also used to cool the vehicle, we assume a 
hydrogen-rich mixture whose stoichiometry is represented by 

2<f> H 2 + 0 2 + 3.76 N 2 -* 2(4>~1)E 2 + 2H 2 0 + 3.76N 2 

where <f> 2 1. The 3.76 nitrogen coefficient is appropriate for air. For 

stoichiometric combustion, with <f> = 1, a simple estimate for the maximum 
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change in stagnation temperature yields 

AT 0 = T 02 “ T 01 S 2000 K (26) 

The maximum heat addition then is 

q = Y 1 ?-! at 0 S 2 x 10 6 J/kg (27) 

where we use air values for *f( = 1.4) and for the gas constant R(= 287 J/kg-K) . 

Table I is taken from reference 12 and provides typical atmospheric data 
at several representative altitudes, h, for hypersonic flight. We shall use a 
nondimens ional heat addition parameter, defined as 

Q - fa (28) 

Based on Table 1 and equation (27), Q = 30 is a representative value for the 
maximum amount of heat addition under stoichiometric conditions. Since T^j is 
given by 

T 01 = T o» = T . (* + ¥ $ ' 29 > 

the stagnation temperature ratio across the combustor can be written as 




y -1 .2 


(30) 


With a given M^ # it is useful to determine the value for Q, denoted as 

Q , for which the flow in a constant cross-sectional area duct will thermally 
K 

choke. This is given by Rayleigh line theory as [see pg. 435 of reference 10 
with Mg = 1] 
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( 31 ) 



T* T 
0 _01 

T m T 
01 00 


r M - j!ii + V "9 

2(t 2 - 1) M^(l + ^ «?) 


t 

where Tq is a reference (at M = 1) stagnation temperature, and may be 
subsonic or supersonic. Figure 5 shows this relation with y = 1.4 for two 
values of M^. Thermal choking is avoided in a constant cross-sectional area 
duct by having Q 0 £ Q. Thus, if M =7 and Q = 30, then M 1 must exceed 5 to 
avoid thermal choking in a duct with a constant cross-sectional area. 

For a large thrust, we are interested in Q values that may appreciably 
exceed Q^. Thermal choking is then avoided by Ag/A^ having a non-unity value 
such that Mg has its prescribed sonic or supersonic value. 

It is also useful to obtain an estimate of the maximum static 
temperature, T^, which typically occurs at the exit of the combustor. This 
estimate is needed for assessing the approximate validity of the non-reacting 
and calorically perfect gas assumptions. For these assumptions to be 
reasonable, T^ should not be too large, otherwise real gas effects become 
very significant. In turn, this limitation restricts M to the lower range of 

GO 

hypersonic flight Mach numbers. This temperature is obtained as 







1 + 1 r 
1 + ¥ ^ 



(32) 


where = T^. With specified values for y and M^, this relation yields an 
upper bound on Q for a given value of T^/T^. 

From references 13 and 14, we see that dissociation of equilibrium air is 
negligible below a temperature of about 2778 K (5000°R) . For air at 30 km 
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(see Table 1), this means that T /T should satisfy 

mx co 

jnx < 2778 _ 19 9 
T m ' 227 

Our parametric study will usually adhere to this limit. On the other hand, 
equilibrium vibrational excitation of the 0 ^ and molecules becomes 

significant at a much lower temperature. For instance, for equilibrium air, y 
= 1.294 at 2778 K (see reference 14). While vibrational excitation is not 

included in the subsequent model, we will make a rough estimate of its effect 
on the thrust. 


Thrust 

As equations (10) and (11) show, the thrust is given by the difference in 
the impulse function at the exit and inlet stations. For the problem at hand, 
the thrust provided by the combustor is 



(33) 


where a caret denotes a dimensional value. For the overall engine, we have 



,F 0 F 




(34) 


For a perfect gas, equation (11) can be written as 

F = pA(l + tM 2 ) (35) 

Since both thrusts are proportional to F^, it is convenient to define a 
normalized thrust as 
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(36) 


r- 


A 

Poo A l 


This normalization is identical to that used in equation (20) for the drag. 
It is important to note that equations (33) and (34) hold even when skin 
friction, a normal shock wave, heat addition, or a gradual change in the 
cross-sectional area are present. 

Impulse Function 

With the aid of 


m = pAV 

T 0 = T ( 1 + Hr 1 « L ) 


2 YP 

3 P 

equation (35) can be written as 


F = 


-(?)* 


L-L.Ig£ 

"l 1 + H 2 rf 


(37a) 


where a is the speed of sound of a perfect gas. As with equation (35), this 
relation is a point function; its use in equations (33) and (34) means the 
intervening flow may contain shock waves, skin friction, an area change, and 
heat transfer. By evaluating the right side of equation (37a) at M = 1 

we obtain 
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(37b) 


F_ 1 + vM 2 

F* J^2(y+1)M^^1 + M 2 ) J 1/2 

♦ 

where we have assumed Tg has the same value for F and F $ 9 Later, this 
assumption is dropped. ^ 

Equation (37b) is shown in figure 6 for y = 1.4. For this y, the 

± * 

asymptotic value of F/F is 1.429 as M 4 «. Thus, F/F has relatively little 
variation when the flow is supersonic. Later, equation (37b) is used for the 
isentropic nozzle flow in the form 


"2 1 +T"! 

M 3 1 + 



(38) 


The double valued behavior evident in figure 6 can be understood by 
writing F as 

F = pA + mV 

where the lass flow rate, m, is a constant. Consider flow in a duct with a 
fixed exit area A. At a low subsonic Mach number, V approaches zero while p 
approaches its stagnation value, thus the pA term dominates. At a high 
supersonic Mach number, the rightmost term dominates, since p approaches zero 
whereas V approaches its limiting value. Note that F approaches infinity as 
M ■* 0. 

The foregoing infinity motivates us to consider an engine with isentropic 
flow that operates as a shock-free diffuser with a supersonic inlet and a 
subsonic exit. When M is large, the engine has a large thrust since a large 

CD 

fraction of the free-stream kinetic energy is converted into internal energy. 
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The positive thrust stems from a low static pressure on the walls upstream of 
the throat, where the flow is supersonic, and a large pressure that 
approaches the stagnation value on the walls downstream of the throat. This 
thrust does not require any heat addition. By the momentum theorem, it is 
unaltered if the coordinate system is fixed with the engine, as in our 
analysis, or the engine is traveling with velocity V into quiescent air. 

The foregoing discussion presumes the flow decelerates isentropically to 
a subsonic value. If the flow remains supersonic everywhere within the 
converging/diverging duct, it is still possible to obtain a (relatively small) 
positive thrust, providing the exit area exceeds that of the inlet. In fact, 
this happens in the subsequent parametric model, i.e., the (inviscid) engine 
has positive thrust even when the heat addition is zero. 

Different definitions for the thrust and drag are common. The 
definitions typically depends on the experimental test configuration or on the 
sophistication of the model. In any case, clear definitions are generally 
required. 


Normal Shock Wave 

Within the context of the analysis, a normal shock wave may be present 
whenever the upstream flow is supersonic. Conditions just upstream and 
downstream of a normal shock are denoted with u and d subscripts, 
respectively. Conservation of momentum across the shock is given by 

(P + pv 2 ) u = (p + pv 2 ) d 

Multiplication by yields 



and the impulse function does not change across a normal shock. 
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In certain circumstances, a change in* location of a normal shock inside a 
duct does not alter the impulse function at the exit plane of the duct. The 

one nontrivial example of this is Rayleigh flow as shown on pp. 123-125 of 
reference 10. 

More generally, however, the location of a shock wave does alter the exit 
plane value of the impulse function, and thereby alters the thrust. This is 
most readily demonstrated for a nozzle flow that contains a normal shock wave 
inside its divergent section. With inlet conditions fixed, the exit plane 
subsonic Mach number depends on the location of the shock. Consequently, the 
value of the exit plane impulse function also varies with the location of the 
shock. 

As evident from figure 6, a large thrust may occur if the exit Mach 
number Mg is very small compared to unity. This will occur if the back 
pressure of the nozzle is sufficiently high so that there is a normal shock 
near station 2 that occurs in either the combustor or the upstream part of 
the nozzle. There will be a high static temperature downstream of the shock 
with thermal dissociation of the normal combustion products, such as HgO, NO^, 
Ng, — . The dissociation process, of course, prevents the static temper- 
ature from escalating to an even higher value. It also represents a 
significant endothermic process that reduces the stagnation temperature and 
the thermal efficiency of the engine. 

In addition to the above discussion, the presence of a normal shock, or, 
more realistically, a shock system, causes boundary- layer separation and high 
local heat transfer rates. In this circumstance, an idealized one-dimensional 
analysis would grossly overestimate the engine’s thrust. 

In view of these factors, the flow inside the engine hereafter is assumed 
to be free of a normal shock wave. 
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Influence Coefficient Method 


Our discussion in this subsection does not presume inviscid flow in the 
combustor, except when stated. In the influence coefficient method, the 
quantities 

dA fb ar dx 
A • T 0 ’ f D 


act as forcing functions, where x is axial distance, D is the hydraulic 
diameter, and is the skin friction coefficient. The above quantities are 
presumed to be known functions of either M or x. Once known, other variables 
are determined by equations such as equations (8.7) in reference 10. Our 
analysis utilizes three such equations: 


dF 

F 


1 dA x M 2 
1 + ym 2 a 2 1 + tM 2 



(39) 


M 2 - 1 dM 2 „ dA 
X+l^M 2 A 



(40) 


M 2 - 1 dg dA 

TM 2 p '~ a 




(41) 


Equation (39) is equation (8.12) in reference 10, while equations (40) and 
(41) are taken from p. 433 of this reference. Although dT^ does not appear in 
equation (39) , the impulse function depends on the heat addition through 
changes in the Mach number. The above relations stem from the one-dimensional 
conservation equations for a steady flow, the state equations of a perfect 
gas, and the definitions of the impulse function and Mach number. 
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If we logarithmically differentiate equation (37a) with dm = 0 and 
substitute the resulting value for dF/F into equation (39), we obtain equation 
(40). Thus, equation (37a) is an exact solution of equations (39) and (40). 
We thus write 


F_ JJl 1 + rM 2 

F l" M i +y m? 


(l + 


1 + 


2 







1/2 


(42) 


for the flow in the combustor. If dA = = 0, we have Rayleigh flow for 
which we can show that = 0* When / 0 and dA = 0, it is the skin 
friction that provides the thrust even though does not appear in equation 
(42). In this circumstance, we can show that ^ irrespective of the 
value of M^. 

We can now write for a combustor with skin friction and a variable cross- 
sectional area 


F 2 Mj i + tm| 


1 1/2 % 
I 

F 1 " ”2 1 + TM 2 

} * V •£ 

T 0! 


1/2 


(43) 


This relation shows the importance of a large Q, or ^02^^01 ’ va l ue f° r 
increasing Fg relative to F^. 

From equation (39), we observe that friction always produces a negative 
dF. The effect of heat addition is not so simple, especially at supersonic 
speeds. While heat addition may increase F because of the stagnation 
temperature ratio in equation (42), it also alters the Mach number. At 
supersonic speeds, the Mach number decreases due to heat addition, thereby 
tending to decrease F. These opposite trends suggest a solution that would 
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maximize the thrust of the combustor. We address this question in the next 
section. 


ANALYSIS 

Thrust 

Within the context of the analysis, the two thrusts of interest are 
written as 




(l + ^ M 2 ) 


1 + I r M V 


Y/(Y-1) 


- 1 


(44a) 






+ m? 


( F _2 F _3 
iF l F 2 



(44b) 


where F^/F^ is given by Eq. (43) and F^/Fg by equation (38). We presume all 
Mach numbers are sonic or supersonic, 'i M^, Mg ? Mg, and Mg may be greater 
or less than M^. 

Several general conclusions can be drawn from equations (44). First, if 


2 

=- > 1 and 
F 1 


M 3 > ”2 


(45) 


we have 

1 ' i _ 3 > 7{_ 2 > ° 

It is possible for ^ ess l^an unity, in which case K 2 is 
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negative. For example, this occurs when 


0 = 0 , 1 $ M 2 < Mj 

An important conclusion is that ^^-2 * s “dependent Mg. (This result 
does not assume an inviscid combustor.) A sonic value for Mg means that a 
relatively large fraction of the thrust is produced by the nozzle. 
Alternatively, if Mg = 
the combustor. 


there is no nozzle and all of the thrust comes from 


In our analysis, various Mach numbers are prescribed rather than area 
ratios. As we have just shown, the overall thrust is independent of Mg and 
therefore is also independent of Ag/A^. Consequently, there is no point to 
maximizing the combustor thrust g, s * nce the overall thrust remains fixed. 

With prescribed area ratios instead of Mach numbers, we see that the 
combustor thrust has a maximum value when Ag/A^ > 1 and the heat addition is 
just sufficient to result in Mg = 1. In this circumstance, the pressure Pg at 
the exit of the combustor has its maximum value. (There is a structural limit 
on the magnitude of the combustor pressure. Because of the high aspect ratio 
of a roughly rectangular combustor cross section and its severe thermal 
environment, this pressure limit is rather low.) In any case, it is the 
overall thrust that should be maximized, not that of the combustor. Shortly, 
we derive a condition for maximizing when the Mach numbers are 
prescribed. 


The independent parameters in this influence coefficient model are y 9 M , 

CD 

Mi, Mg, Mg, and Q. Once these are specified, with the exception of Mg, then 
^2 is determined. Since the combustor is located inside the vehicle, it is 
advisable to keep Ag/A^ from becoming unduly large. As will be shown, this is 
achieved by setting Mg = 1. (We also discuss cases when Mg > 1 and Ag/Aj is 


25 



not unduly large. ) 


It is useful to write equation (44b) as 


£3 = (‘ + 


(i ♦ ** ifl 


2 ”00 


1 + V ”ij 


T/(T-1) 


Mj 1 + ym| 


1/2 „ 

M 3 1 + yMj 

+ ¥ 4 

T 01 


1/2 


- 1 


(46) 


We suppose M^, Mg, and are kept fixed, and examine the dependence of 

Yi- 2 on Mj. We, therefore, write this relation as 


T 

i-3 


BM, 




. (T+1)/[2(y-1)] 

“lj 


c(i + tmJ) 


(47) 


where the fixed parameters B and C are 


B = 





1/2 


c * p ¥ « 2 J t/<t " ) 

An extremum value of c ? is found by differentiating equation (47) with 
respect to M^, with the result 
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K- l ) 

TCM, 


1(1 * ¥ «r 2 J 


By setting the right side equal to zero, we obtain 


M 1 = 1 


which yields a maximum value, CT . given by 

i"o ,mx 

2 ><t+«/[ 2 (t-i)J, «r/(H) 


X 


1-3, mx Vf+1 


ft) (* + V 


J ,( 1 + r $ 


M 1 + ¥ < 



p2\ 

j 1/2 

* T or 


(¥ /2 “ 


( 48 ) 


(49) 


The other root 

( J + *T "! ) VZ B = tcMj 

yields a negative value for K g, and therefore is of no interest. 

From equation (46) , we observe that the largest permissible values for 
M^, Mg, and Tgg/Tgj yields a maximum value for c^^ 3 * The resu lt f° r ^3 i s 
evident from equation (42) and figure 6 . These parameters, however, are 

constrained by a number of conditions, such as a maximum combustor static 
temperature. 

Area and Pressure Ratios 

We obtain the pressure and area variation under the assumption of 
inviscid flow. For the combustor, this is accomplished by integrating 
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equations (40) and (41) with = 0. The results for the nozzle integration, 

where we also have dT^ = 0, are provided by the isentropic results, equations 

(24) and (25). Aside from isentropic flow, there are two other simple special 

cases for the combustor. In the first of these, we assume dA = 0. This 

yields the Rayleigh flow solution, which requires Q £ Q. Neither this 

H 

solution nor the isentropic one are of interest for the combustor. 

The last special combustor case, dM 2 = 0, easily yields for the area 

ratio 




dnH 2 )/2 


(50a) 


With the aid of equation (41), the pressure variation is 


P 1 T 01 


(50b) 


For Y = 1*4 and M = 1, the area grows slightly faster than linear with Tq/Tq^* 
For moderate supersonic Mach numbers, the rate of growth can still be 
relatively slow when t q/ t qj^ does not rapidly increase with x. 

We perform the integration of equation (40) for the combustor by 
arbitrarily assuming a relation between and M. The equation to be utilized 
is 



(51) 


This relation is chosen for its analytical simplicity. Other relations can be 
used, e.g. , one for A(x) and another for T^(x). It is important to note that 
the combustor and overall thrust values are unaltered by any of these choices. 
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Equation (51) yields 


5 

T 0 


= OtdM 


( 52 ) 


where 


« = 


1 Mo T<!2 


II || ’* il rp 

M 2 M 1 *01 


(53) 


and T Q2 /T 0 i is provided by equation (30). With equation (52) and = 0, 

equation (40) becomes 


M 2 - 1 cM 2 


= 2 r - «(i + ym 2 )<w 


l+l^M 2 A 


The leftmost term integrates to 


f ^ M 2 - 1 dM 2 _ jj 

/VI 


I+l-l 

Y-l 


Im j | 




We thereby obtain 


A_ 

A, 


5 

M 


1 + 3C-1 m 2 ^ 


1 + 3C r< 


I+L 


2(y-D 


exp {l [(” ~ M l) + 3 K " ”l)]} 


(54) 


for the variation of the combustor’s cross-sectional area with Mach number. 
[The variation of A with axial distance requires specification of M(x) or its 
equivalent. J 

The pressure ratio is evaluated by eliminating the dA terms in equations 
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(40) and (41) , with the result 


dg - _ 1 Im 2 ^ dM 2 ) 

P ' " 2 I T 0 l^H 2 


With equation (52), this becomes 


? = | KM2dM + 


dM 2 




which integrates to 


E_ 

p i 




1 


r/(r-i) 

exp I- 


( 55 ) 


(56) 


Observe that equations (54) and (56) reduce to isentropic relations when 
« = 0. We attain Ag/A^ and p^/p^ by simply setting M = Mg in these equations. 
In this instance, equation (54) reduces to 



1/2 + t(Mi + M 1 M 2 + m|)/6 


(57) 


with the aid of equation (53). Here, the is subscript denotes an isentropic 
point relation. Note that equation (57) reduces to equation (50a) when M ^ = 
Mj. Actually, equations (51)-(54) and (56) are indeterminate when in 

which case they are replaced by equations (50). 

A useful result is obtained by considering and as fixed, with 

*5 
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the additional restrictions: 


M 1 Z 1 , Mg 2 1 



1 


In this situation, (Ag/A*) ^ and the (^Qg/T^) ^ ac ^or in equation (57) 
increase with Mg. Thus, 


V A i 


is a minimum when 


«2 = !• 


with the value 


tz 

A, 


Tq 2 %l/2 + T(Mf + M 1 + l)/6 



(58) 


While both numerator and denominator increase with M-^ in this 
numerator increases more rapidly than does the denominator, 
exceeds unity. 

A Brayton cycle is obtained by replacing equation (52) 
Equation (55) then integrates to 


relation, the 
Hence, Ag/A^ 

with dp = 0. 


f T ° fo . _ f ** dH 2 

T 0 1 T ° M? ^ 


!o_ 1 + *r f 

T 01 M 2 1 + 


The heat addition Q can no longer be prescribed, since Tgg/Tgj is determined 
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. Alternatively, if Q is prescribed, then Mg cannot be 
prescribed. In this circumstance, thermal choking must be avoided. 

In this analysis, we prefer the greater flexibility of prescribing both 
Mg and Q. Our results frequently correspond to a slightly favorable pressure 
gradient for the combustor, i.e., 1° ^^ ese cases * a Brayton cycle 

would result in a reduced value for the Ag/A^ area ratio and frequently would 
result in thermal choking. 

PARAMETRIC RESULTS 
Nominal Cases 

After several preliminary investigations, three nominal cases were 
chosen. These cases have Y = 1.4 and the values shown in table II. The 

nominal cases are designed to evaluate the effect of changing M with the 

GO 

pressure ratio Pg/P^ set equal to unity. This is accomplished by using 

Y/(T-1) 


£3 £2 £l 

p 2 P 1 p a. 


1 + n 
Z qo 


1 + 




-TM2/2 


02 

r 0L 


= 1 


by setting M = Mg 


where equation (50b) provides Pg/Pj- This relation then yields 


«3 = 




-(Y-1)m|/2 



1/2 


(60) 


Table III lists results for these cases. We see that Q 0 exceeds Q, 

except for the third case, and all cases adhere to an upper bound of 12.2 for 

T^/T^. With increasing M^, this upper bound can be met only by increasing 

and decreasing Q. Thus, at = 9, ^Q2^01 * s on ^y 1*166 although Q is still 

1/3 of its value at M =5. In case 3, combustor operation is either very 

00 
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fuel or oxidizer rich. 


Although Mj increases with decreases because of the decrease in 

Tog/Tgi* The values shown for A^/A^ are rather modest and should be easily 
achieved with a practical combustor. The rapid increase in A^/Ag is expected, 
in view of the increase in M^. All three combustors have a modest favorable 
pressure gradient, which should help reduce boundary-layer separation 

_3 

phenomena. The values in the p^/pg column are to be read as 6.66 x 10 ; they 

are as expected. 

The final two columns show the combustor and overall thrusts. All values 
greatly exceed the earlier drag estimate of 206. The fraction of the overall 
thrust provided by the case 1 combustor is 52.5%. This large fraction occurs 

it is a result of the higher pressure 
level in the combustor relative to the nozzle. The fraction falls rapidly 
with because of the decrease in A^/Aj. Hence, for case 3, the fraction is 
only 23.4%. 

The rapid increase in both thrusts with M occurs despite the decrease in 

00 

Q. The increase is due to the isentropic compression upstream of the 
combustor. This is evident from the multiplicative factor 


even though A^/A^ greatly exceeds A^/ 


1 + 


1 + 



(61) 


that appears in equations (44). This compression corresponds to a significant 
value for the cross-sectional area ratio A /A,. For cases 1 and 3, this ratio 
is 25 and 194, respectively. 

As is evident from table II, we have used a constant Mach number process 
for the combustor. This does not necessarily represent an optimum condition, 
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it merely seems to provide a better compromise for meeting the various 
constraints, such as a reasonable value for Ag/A^ and a favorable pressure 
gradient, than alternative choices. This combustor process, however, does 
provide a thrust optimum when = 1. 

For a later comparison, we need the thrust and mass flow rate m 

scaled by the exit area of the nozzle. These parameters are given by 
A 



h ^2 

^00 Ag Ag ^1-3 


(62) 


m 

A„ 


A, An 

p a Mi r r 
00 * 1 h. H 


fl + 3C-1 M 2 ) 


1 + I r M i) 


(T+1)/[2(T-1)] 


(63) 


where the free stream speed of sound, a , and p are provided by table I. We 

00 GO 

use the nozzle’s exit area for the scaling, since it may be regarded as a 
measure of the vehicle’s cross-sectional area. With M =7 (case 2) and an 

CO 

altitude of 30 km, we obtain 

A 



= 2.58 x 


10 4 Pa 


T- = 25.7 kg/m 2 -s 
A 3 


(64) 


Since one atmosphere pressure at sea level is 1.015 x 10 Pa, the thrust level 
appears to be rather low. This value is placed in better perspective in the 
final section. 


, and Mg 

Table IV shows the effect of varying Mj, Mg, and Mg. All cases use y = 

1.4 and the nominal cases appear as cases 1.0, 2.0, and 3.0. Variations about 

case 1.0 are listed as l.xx, where xx = 1,2,..., and these cases have M =5 

00 


Influence of Mp 
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and Q = 30. A similar statement holds for cases 2.xx and 3.xx. With the 

exception of cases 2.3 and 3.3, all the cases have a T /T below 12.2. The 

mx oo 

two exceptions exceed the T /T limit because M, = 1. Most cases have M, = 

mx oo i 1 

MgJ the few exceptions are sufficient to illustrate the trends that occur when 
this condition doesn’t hold. 

With the exception of cases 1.2 and 1.3, all cases have modest values for 
Ag/Ai* The two exceptions demonstrate that Ag/A^ * s sens i^ive to Mg when M^ 
is relatively small. While there is a wide variation in the Ag/Ag values, 
these values are nevertheless expected. 

With one exception, case 1.2, all Pg/p^ va l ues range from 0.187 to 0.9 
and a favorable pressure gradient is present in the combustor. As with Ag/Ag, 
Pg/Pg has its expected values. Only the nominal cases have Pg/p ^ = 1. A 

number of cases have a p 0 /p value that is below unity. Only in case 1.2, 

u oo 

however, is the value so low that boundary- layer separation inside the nozzle 
would be anticipated. 

The last two columns show that the thrust is not overly sensitive to 

small changes in when M is small. When M = 9, however, there is 

1 00 00 

considerable sensitivity. As shown by cases 1.3, 1.4, 2.5, and 2.6, this 
sensitivity is associated with not Mg. As previously noted, an M^ value of 
unity yields a maximum value for while a change in Mg alters 

Ag/Ag, and the fraction of the thrust produced by the combustor but not the 
overal 1 thrust . 

Examination of cases 1.5, 1.6, 2.7, 2.8, 3.5, and 3.6 show that modest 

changes in Mg, as expected, have little effect on the thrust. It is worth 

mentioning that a modest decrease in Mg from its nominal value produces a 
significant decrease in Ag/Ag. Thus, a considerable truncation of the nozzle 
appears to be warranted. 
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The Mj and Mg values in table IV are relatively small and don’t differ 
appreciably. It is of interest, therefore, to consider a larger Mj value and 
a wider range of Mg values. This is done in table V, where we fix 

T = 1.4 , M =9 , = 4 , M_ = 8 

GO I O 

The first three cases have Q = 10, which is the nominal value for this M 

00 

value, while the last three cases have Q = 20. Within each Q group, T /T , 

* mx co 

A 3 /A 2 , and Pg/p 

i ncrease rapidly with Mg. Observe that is a constant within each 

group. Cases 3.7 and 3.10 have an adverse pressure gradient in the combustor 
even though Ag/A^ exceeds unity. As noted, Vx- 2 increases rapidly with Mg to 
the point where most of the thrust is generated by the combustor. In this 
circumstance; however, the large Ag/A^ value may not be practical. It is 
important to note that the overall thrust level in this table is an order of 
magnitude less than that in table IV. This is a consequence of the larger 
value. 


j fall rapidly with Mg. On the other hand, 


V 4 r 


p 3 /p 2 . 


and 


Influence of Q and Y 

The effect of the final two parameters is shown in table IV. Except for 
Q and Y> all l.xx cases have the prescribed values of the 1.0 nominal case. 
Examination of the T /T column shows a considerable Q effect when M is 

mx ai oo 

small and a much smaller effect when M is large. The reason for this is that 

co 

most of the temperature increase, when M is large, is due to the upstream 

GO 

isentropic compression. 

The changes in 3 with Q are as expected. One interesting observation 
is the substantial thrust that occurs when Q = 0. In this case, the flow is 
entirely isentropic as noted earlier. Undoubtedly, a more realistic model 
that incorporates shock waves, skin friction, etc. would alter this finding. 
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Our final observation is the sensitivity to even a modest decrease in Y, 

especially when M is large. The increase in the thrust is entirely due to 
00 

the compression upstream of the combustor as represented by the item (61) 
factor in equations (44). This result indicates the importance of evaluating 
real gas effects in the forebody and scramjet inlet regions. These effects 
not only include compositional changes but also vibrational excitation of air, 
which is significant at relatively low temperatures. 

SUMMARY AND DISCUSSION 

We first compare the thrust and mass flow rate values of the preceding 
section with that of a conventional rocket engine. For a meaningful 
comparison, this is done on the basis of the nozzle’s exit area. For a 
rocket, we use the V-2 engine, whose performance has certainly been surpassed. 
However, data (see reference 15) for this engine is readily available, and 
results in 

, ^ = 274 kg/m 2 -s 

where A^ represents the exit nozzle area of the V~2 engine. These values are 

contrasted with the scramjet values in equations (64). We observe that the V- 

2 thrust and mass flow rate are 21 and 11 times that of the scramjet, 

respectively. Of course, the V-2 engine must overcome gravity for the 

missile, as well as a small missile drag. On the other hand, the scramjet 

must simply overcome drag and for this it suffices. Nevertheless, it is 

useful to note that the scramjet is a low power density engine, as measured by 
A 

^l-3 /A 3‘ 

Another commonly used scramjet performance parameter is the specific 
impulse, more properly referred to, as it sometimes is, as specific thrust. 
It is defined as 


r 5 

-Z- = 5.46 x 10° Pa 
a 3 
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z 

■f 


where m ^ is the mass flow rate of the fuel. This later parameter has not 
appeared in the analysis because we have not evaluated m_p. 

The first three sections of this report are primarily tutorial. They 
establish the necessary background for a simple quasi-one-dimens ional model 
for estimating scramjet performance, which is the subject of the preceding two 
sections. Among the various assumptions that are utilized, two standout as 
limiting the scope of the model. These are the assumptions of isentropic flow 
upstream of the combustor and of a calorically perfect, non-reacting gas. 
These assumptions tend to limit the validity of the model to relatively low 
hypersonic flight Mach numbers. 

Within the constraints of the model, a number of conclusions can be 

drawn. These may be summarized as follows: 

1. The overall thrust is a maximum when M , M OJ and Q are maximum and 

00 d 

when M, - 1. Of these conditions, a large M value and M, = 1 are the most 

1 co x 

important. When Mg is already large compared to unity, any further increase 

in Mg provides only a minimal thrust increase. While Q is more important than 

M 0 , its effect is less important than M, or M . As M is increased, starting 

from a value of about 5, M^ must be increased above unity if the upper limit 

on T /T is to be maintained. In this circumstance, a constant Mach number 
mx oo 

combustion process is found to be satisfactory, although the total thrust 
falls rapidly with increasing M^. 

In view of the thrust’s modest dependence on Q, we believe real gas 
phenomena, such as vibrational excitation and molecular dissociation, needs 
further study. It also may be advisable to consider T /T values greater 

mx oo & 
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than 12.2, especially at large values of M^. 

2. The fraction of the thrust associated with the combustor is 
substantial, often exceeding 50% of the total thrust. This fraction depends 
on the value of Mg- A value of Mg = 1 provides a minimum value for the 
combustor area ratio Ag/A^ and a maximum for the fraction of thrust produced 
by the nozzle. 

With Mj and Mg prescribed, it is important to note that the total thrust 
is independent of Mg, even when the combustor flow is viscous. For an 
inviscid combustor, the thrust is independent of the heat addition 

distribution inside the combustor. 

3. In view of items 1 and 2, the M = 1 relation for the combustor 

represents an optimum thrust condition. In this circumstance, the thrust is a 

maximum and the combustor area ratio is a minimum. This condition represents 

an area ruling relation for the combustor. In other words, when the variation 

in stagnation enthalpy is known as a function of axial distance, the 

combustor’s cross-sectional area, A(x), can be determined. As noted, when M^ 

increases, the value of T /T becomes a constraint that requires M., to exceed 

mx ao 1 

unity. In this circumstance, the combustor condition, dM = 0, appears to be a 
suitable alternative that provides a modest favorable pressure gradient for 
the combustor. 

4. The thrust is sensitive to the value of the ratio of specific heats. 
This effect is due to the large compression that occurs in the forebody and 
scramjet inlet. This result emphasizes the importance of properly assessing 
real gas phenomena in the flow upstream of the combustor. 
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TABLE I. Standard Atmospheric Data 12 


h, altitude 

Poo f 
Pa 

Poo 9 

kg/m 3 

Too 5 

K 

CO * 

m/s 

km 

kft 

25 

82.0 

2.549 x 10 3 

4.008 x 10 -2 

221.5 

298.4 

30 

98.4 

1.197 x 10 3 

1.841 x 10" 2 

226.5 

301.7 

35 

115. 

5.746 x 10 2 

8.463 x 10“ 3 

236.5 

308.3 


TABLE II. 

Nominal 

Cases 

with y = 

1 . 4 and 

(P 3 /Pj = 1 

Case 

M 

oo 

M i 

M 2 

M 3 

Q 

1.0 

5 

i 

1 

4.486 

30 

2.0 

7 

1.5 

1.5 

6.293 

20 

3.0 

9 

2 

2 

8.429 

10 
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